Stability of magnetic vortex in soft magnetic nano-sized circular cylinder 
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Stability of magnetic vortex with respect to displacement of its center in a nano-scale circular 
cylinder made of soft ferromagnetic material is studied theoretically. The mode of vortex displace- 
ment producing no magnetic charges on the cylinder side is proposed and the corresponding absolute 
single-domain radius of the cylinder is calculated as a function of its thickness and the exchange 
length of the material. In cylinders with the radii less than the single-domain radius the vortex 
state is unstable and is absolutely prohibited (except if pinned by material imperfections), so that 
the distribution of the magnetization vector in such cylinders in no applied magnetic field is uniform 
(or quasi-uniform). The phase diagram of nano-scale cylinders including the stability line and the 
metastability region obtained here is presented. 



The small magnetic nano-scale cylinders made of soft 
magnetic materials recently gained attention due to the 
progress in fabrication and observation techniques and 
also because of their possible applications in magnetic 
random access memory (MRAM) devices. Within such a 
cylinders of circular shape in a certain range of sizes the 
magnetic vortices are frequently observed (see e.g. 
^, Q] ) . For applications, which usually try to avoid vortex 
formation (such as MRAM cells, Q), it is important to 
know the sizes of the cylinder where the vortices do not 
form. 

In thin ferromagnetic cylinders with the thickness L 
of the order of a few Le (where Le = \JC/Mg is the 
exchange length, C is the exchange constant of the mate- 
rial, Ms is the saturation magnetization) distribution of 
the magnetization vector can be assumed uniform along 
the cylinder axis. Then, there are two characteristic sizes 
of the cylinder important for the presence of the vortex 
state versus the uniformly magnetized one in zero applied 
magnetic field. The first is the single-domain radius Req, 
which is defined as a radius of the cylinder (at a given 
thickness) in which the energies of the uniformly magne- 
tized state and the state with the vortex are the same 
1^. In cylinders with radii R < Req the uniformly mag- 
netized state has a lower energy than the vortex state. 
However, the metastable vortices still may be present in 
cylinders with radii below Req- There is another char- 
acteristic radius, the absolute single domain radius Rs of 
the cylinder, which is obtained from the requirement that 
the vortex is unstable (therefore is absolutely prohibited) 
in cylinders with R < Rs- 

The rigorous calculation of Rs requires evaluating the 
second variation of the energy functional including the 
long-range dipolar interactions, which is currently be- 
yond possibilities of analytical methods. The other way 
to estimate the stability radius is to assume the precise 
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way (mode) the vortex loses its stability and then to cal- 
culate the stability radius with respect to that process. 
There is an infinite set of possible candidate modes. Pro- 
vided calculation of the energies is rigorous, the result for 
a particular candidate is a lower bound for the stability 
radius. The "true" Rs is the highest of all the lower 
bounds for all possible modes. This is also easy to un- 
derstand from a thought experiment of having a cylinder 
with vortex and slowly shrinking its radius, as soon as 
the largest absolute single domain radius is reached the 
vortex loses stability according to that particular mode. 

The first estimate of the absolute single-domain radius 
Rs was obtained for the circular cylinder by Usov and 
coworkers Q for the mode corresponding to the uniform 
translation of the vortex. Such a translation is accom- 
panied by formation of magnetic charges on the cylinder 
sides. 

In small and flat ferromagnetic cylinders of isotropic 
(ideally soft) material the following hierarchy of energies 
is present (by decreasing their importance for mini- 
mization): exchange energy, energy of magnetic charges 
of the cylinder faces (because the cylinder is fiat), energy 
of magnetic charges on cylinder sides, energy of volume 
magnetic charges. The energy of volume charges is the 
least important one, because of a general amplification 
of the effects associated with surface compared to the 
ones associated with volume, as the size of the parti- 
cle decreases. Thus, it may be expected (and further 
calculations in this work support this conjecture) that 
a mode with non-uniform vortex deformation having no 
side magnetic charges (but having volume ones) may be 
more favorable for the vortex stability loss. 

To find such a mode is not an easy task because once 
the vortex is deformed the problem loses the cylindri- 
cal symmetry and becomes two-dimensional (leading to 
a system of partial differential integral equations) instead 
of one-dimensional (described by a single ordinary differ- 
ential equation) if the axial symmetry is present. It is 
possible to write the solution for the displaced magnetic 
vortex with no side magnetic charges and very small face 
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charges at the vortex center exactly minimizinff the ex- 
change energy functional in the following form [p| 
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where z = X + lY, i = y^—l, X and Y are Cartesian 
coordinates in the cylinder plane, so that the axis Z is 
parallel to the cylinder axis, the complex parameter a 
describes the vortex displacement (for a — Q the vor- 
tex is centered), the real parameter c is related to the 
vortex radius, bar over a variable denotes the complex 
conjugation, the components of the magnetization vec- 
tor M{X, Y) are given by Mx + iMy = Ms2w/(1 + ww) 
and Mz — Ms{l — ww) / ( 1 + ww) . An example of config- 
uration (|l|) is shown in Fig. |^. In the case of no applied 
field (the case when the field is non-zero was considered 
elsewhere the phase of the complex number a is not 
important, this parameter will be considered real in the 
rest of this work. The expression (|l|) is not arbitrary, it 
follows from the analysis performed in Q, that it is the 
only way to displace the vortex center so that it's struc- 
ture keeps minimizing the exchange energy functional ex- 
actly and has no surface magnetic charges on the cylin- 
der side. The magnetization distribution (^ has surface 
magnetic charges on the faces of the cylinder localized 
near the vortex center, which is a topological singularity 
Unlike side charges, the charges on the faces can not 
be avoided if the cylinder is in a state of non-uniform 
magnetization (topologically charged state), such as vor- 
tex. 

The values of the parameters a and c need to be found 
from minimization of the magnetostatic energy for the 
configuration (^). It is clear that in absence of magnetic 
field, the equilibrium a = due to symmetry, in this 
case (^ reduces to the structure of an axially symmetric 
centered vortex exactly coinciding with the one studied 
by Usov|^. The parameter c = Ry/R, where Ry is the 
equilibrium vortex radius calculated in ^ . 

The stability radius can be found from the condition 
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which means that the energy minimum on the variable 
a turns to maximum. In this expression E^x and Ems 
are the total exchange and magnetostatic energy of the 
cylinder respectively. 

The normalized exchange energy of the particle with 



magnetic configuration (^ is 
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where V — LttR^ is the particle volume. 

Due to the lack of space the change (with a) of the 
vortex core magnetostatic energy (face charges) will be 
1 
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FIG. 1: The distribution of the magnetization vector (^) 
corresponding to a displaced vortex with c = 1/2, a = 1. 
Only projection of the magnetization to the cylinder face is 
shown, shorter arrows mean non-zero out-of-plane magneti- 
zation component. 



neglected in this paper. This approximation is equivalent 
to neglecting terms of the order c'^ <C 1. The density of 
volume magnetic charges (also neglecting the fact that 
inside the vortex core |/(z)| < 1 they are reduced, equiv- 
alent to throwing out terms of the order ^ 1) is 
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in the polar coordinate system r, centered at the par- 
ticle axis. Their normalized energy is 
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where straightforwardly 
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— . - equilibrium radius of vortex core R^/L 




equilibrium single-domain radius R| 
absolute single-domain radius 
previous estimation of R^ 
asymptotic of the R^ at L/Lg«l 
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which can be simplified into 
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FIG. 2; The phase diagram of nano-scale cylinders. The re- 
man numbers mark the regions separated by the solid lines 
corresponding to different ground states of the particle ^, |^ : 
I - vortex state, II - magnetized uniformly in-plane, III - mag- 
netized uniformly parallel to the cylinder axis. The shaded 
region marks the range of cylinder's geometrical parameters 
where the vortices are metastable, below that region no com- 
plete vortices can be present. The previous estimation of Rs 
corresponds to the uniform vortex displacement Q 
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Equation (0) was then solved numerically for Rs/Le at 
different L/Le, the result is shown in Fig. It is possible 
to obtain the asymptotic of this dependence at L <C Le 
analytically by expanding (x << 1) = kx"^ + 0{x^), 
where k = 0.0882686 ... we get 
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This asymptotic was used to verify the results of the nu- 
merical computation and is also plotted in Fig. |[ 

Neglecting the vortex core is not such a bad approxi- 
mation as it seems, in fact, it only breaks when c = Ry / R 
is close to one, which happens near the intersection of 
the Rs and Ry lines in Fig. ^. Including the core into 
consideration is straightforward but involves more com- 
plicated formulas. Such more detailed calculation will be 
published in a forthcoming paper. 
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